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Дмитрий Речкин 

Невидимая динамика солнечной активности 


1. Материал и методы 

1.1. Наблюдаемые циклы солнечной активности, числа Вольфа 

Теперь уже ни для кого не секрет, что центральное тело нашей Солнечной системы, звезда по 
имени Солнце 1 , живет своей собственной, не всегда понятной для нас, но очень активной и на¬ 
сыщенной жизнью. Желающие изучать эту активность неизбежно столкнутся с огромными объ¬ 
емами накопленных данных (прежде всего - многовековыми данными астрономических на¬ 
блюдений). Рудольф Вольф 2 , заинтересовавшийся открытой Генрихом Швабе 3 периодичностью 
количества наблюдаемых солнечных пятен, собрал все доступные с 1610 года данные, и до¬ 
полнил их собственными наблюдениями. В 1848 году он предложил собственный метод расче¬ 
та числового показателя солнечной активности; получаемые по этой методике значения полу¬ 
чили название чисел Вольфа. В 1849 году официальную публикацию чисел Вольфа начала Цю¬ 
рихская обсерватория. 
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Рис.1. График среднемесячных чисел Вольфа с 1750 года 

Р.Вольф первым вычислил периодичность цикла солнечной активности, определив его в 11,1 
лет. В дальнейшем стало понятно, что цикл не всегда является строго «одиннадцатилетним», 
так, в XX веке продолжительность солнечного цикла оказывается ближе к 10,5 годам. Совре- 
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1 В.Р.Цой «Звезда по имени Солнце» - М: «Эксмо-пресс», 2002 

2 Вольф, Иоганн Рудольф (нем. ]оМапп КисІоН" ѴѴоІГ), 7 июля 1816 - 6 декабря 1893; швейцарский астроном и ма¬ 
тематик; профессор астрономии Бернского университета (1844), директор Бернской обсерватории (1847), про¬ 
фессор астрономии в Университете и Высшей технической школе Цюриха (1855); член Французской академии 
наук (1885). 

3 Швабе, Самуэль Генрих (нем. ЗатиеІ НеіпгісГі ЗсИѵѵаЬе), 25 октября 1789 - 11 апреля 1875; немецкий аптекарь, 
ботаник и астроном; систематически вел наблюдения за солнечными пятнами, в 1843 году обнаружил законо¬ 
мерность в изменении их числа; открытие получило широкую известность только благодаря Александру Гум¬ 
больдту, включившему данные наблюдений Швабе в третий том своего труда «Космос» (1851); лауреат Золотой 
медали Лондонского королевского общества (1857), иностранный член Лондонского королевского общества 
(1868). 







менные исследователи говорят не об «11-летнем», но скорее о «22-летнем» цикле активно¬ 
сти, напоминая, что после каждого 11-летнего периода активности магнитные полюса Солнца 
меняются местами, и к «исходному состоянию» наше светило возвращается не каждые 11 лет, 
а примерно каждые 22 года. Поэтому в расчетах используется концепция 22-летнего цикла ак¬ 
тивности Солнца. 

Исследование динамику солнечной активности (пусть и по косвенным, усредненным показате¬ 
лям) математическими методами ставит цель наглядно продемонстрировать возможности та¬ 
ких методов, а также поставить вопросы о применимости метода к исследованиям динамики 
других показателей, полученных прямыми астрономическими, астрофизическими и радиоас¬ 
трономическими наблюдениями (измерениями). 

1.2. 5ТАТІС/ЕХ - калькулятор многомерной статистики 

Идея многомерных статистических расчетов на основе матрицы данных, образованной строка- 
ми-объектами и столбцами-признаками является настолько общепринятой, что не нуждается 
ни в установлении приоритета, ни в каких-либо формальных доказательствах. Необходимость 
построения сложных алгоритмов обработки данных последовательным применением извест¬ 
ного набора стандартных преобразований привела к разработке языка ЯКОБИ 4 , в основе кото¬ 
рого лежат принцип директив выбора преобразований, применяемых к матрице данных, а 
также оригинальный способ указания выборки объектов и/или признаков (фрагмента матрицы 
данных), к которым применяются определяемые директивами преобразования. 

Директивы ЯКОБИ представляют собой фразы на языке, приближенном к естественному, что 
делает понятной схему формирования и обработки данных. Лаконичный способ вывода ре¬ 
зультатов вычислений позволяет представить их как в численном (в виде таблиц и матриц), так 
и в графическом (псевдографическом) виде, что максимально облегчает их восприятие и даль¬ 
нейшее использование. 

Первая реализация ЯКОБИ для электронно-вычислительной машины (ЭВМ) БЭСМ-6 была вы¬ 
полнена В.М.Ефимовым. Необходимость расширения круга пользователей привела к реализа¬ 
ции ЯКОБИ для мини-ЭВМ «Электроника-60» (клон РОР-11/05, Оіёііаі Ециіртепі: Сотрапу). Па¬ 
кет 5ТАТІС 5 обеспечивал необходимые вычисления, несмотря на микроскопический по совре¬ 
менным меркам объем оперативной памяти (56 килобайт!) и практически полное отсутствие 
вспомогательных средств (таких, например, как средства работы с электронными таблицами, 
обычные для офисных систем). 

В настоящее время реализация 5ТАТІС/ЕХ является общедоступной 6 7 , она выполнена в среде 
іагагиз 1 для Бгее РазсаІ Сотрііег (РРС 8 ) и опубликована 9 ' 10 в виде полного комплекта исходных 
текстов. Нумерация версий - сквозная с момента первой реализации 5ТАТІС, в настоящее вре¬ 
мя актуальна версия 1.09, которая и применялась для выполнения всех расчетов, положенных в 


4 В.М.Ефимов, Д.В.Речкин «ЯКОБИ - входной язык пакетов прикладных программ статистической обработки био¬ 
логических данных» - науч.-техн.бюл./ВАСХНИЛ, Сиб.отд-ние, вып.48, с.12-17. - Новосибирск: 1985 

5 Д.В.Речкин «Реализация входного языка ЯКОБИ для мини-ЭВМ "Электроника-60". Пакет 5ТАТІС» - науч.- 
техн.бюл./ВАСХНИЛ, Сиб.отд-ние, вып.48, с.18-24. - Новосибирск: 1985 

6 Текущая реализация 5ТАТІС/ЕХ подробно описана здесь: ІЩрѴ/егоирз.іо/а/зІаІіс 

7 Багагиз: ТНе ргоРеззіопаІ Ргее РазсаІ РАО ЮЕ - НЦрУАл/ѵлл/ѵ.Іагагиз-ісІе.оге 

8 Ргее РазсаІ: Ореп зоигсе сотрііег ІЪг РазсаІ апсі ОЬ|ес1 РазсаІ - НЦрУАлл/ѵѵл/ігееразсаІ.оге 

9 ИнрУ/зоигсеЩгее.пеі/ргоіесіз/іасоЬі-зіабс 

10 М:1:р5://и.1:о/ѵТ3500 








основу данной работы. Калькулятор 5ТАТІС/ЕХ является консольным мульти-платформенным 
приложением (здесь использована сборка для платформы УѴіпсІоѵѵз х86-32Ьіі). 


2. Традиция: периодические колебания и скользящее среднее 

Как известно (см. Рис.1), за период с 1700 по 2019 год наблюдалось 24,5 цикла солнечной ак¬ 
тивности. Иными словами, если исходить из этих данных, продолжительность цикла составит 
13,06 года, что, вообще говоря, сильно отличается от ожидаемого значения. Поэтому для рас¬ 
чета интересующей нас продолжительности цикла применим другой подход. 

Рассчитаем корреляции между текущим годом солнечного цикла и предшествующими годами; 
для того, чтобы свести к минимуму погрешность вычислений, будем отслеживать корреляции 
на протяжении как минимум 150 лет наблюдений (не менее 11-12 наблюдаемых циклов). 

Расчет выполнялся по программе и/о//-/олд./ас (см. Приложение П2.1). Результаты показаны в 

Табл.1: 


Табл.1. Влияние предшествующих лет в циклах солнечной активности 
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Интерес представляют только положительные корреляции, поскольку они помогают опреде¬ 
лить соответствие максимума активности в текущем цикле с максимумами в предыдущих цик¬ 
лах. Для облегчения восприятия такие коэффициенты корреляции выделены жирным шриф¬ 
том. 

Очевидно, вычисленные корреляции показывают (за вычетом 150 начальных лет) наличие 13 
циклов, с продолжительностью (в годах): 11,11,10,11,11,11,11,12,11,10,11,11 и 11 лет (все¬ 
го 142). Таким образом, средняя продолжительность цикла оказывается равной 142/13 = 10,92 
года. 

Традиционные представления о динамике солнечной активности базируются на том, что: (а) ак¬ 
тивность является периодической (в соответствии с тригонометрической функцией типа 5/л или 
Со$), (б) «истинная» активность представляет собой огибающую, сглаживающую месячные (го¬ 
дичные) различия соответствующих чисел Вольфа. 

Аргумент тригонометрической функции за период, соответствующий циклу активности, изме¬ 
няется от 0 до 2п, соответственно, годичное приращение аргумента составит для традиционной 
модели 2тх/10,92 = 0,5754 с достаточной для расчетов точностью. Скользящее среднее вычис¬ 
ляем за 10 предшествующих лет и текущий год (округляя полученную среднюю продолжитель¬ 
ность цикла 10,92 до 11). Полученное значение скользящего среднего (для каждого года) мо- 













дулируем значением тригонометрической функции, ожидая получить картину, похожую на по¬ 
ведение значений чисел Вольфа за те же годы. 

Расчет выполнялся по программе \л/оІ/-ІгасІ.іас (см. Приложение П2.2). Результаты в графиче¬ 
ской форме показаны на Рис.2: 
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Рис.2. Исходный ряд среднегодовых чисел Вольфа и традиционная модель 


Черная линия показывает динамику солнечной активности (среднегодовые числа Вольфа) с 
1700 по 2019 годы, оранжевая - отображает график скользящего среднего (с учетом 11 лет, те¬ 
кущего и 10 предшествующих). Синим пунктиром показан результат моделирования динамики 
модуляцией скользящего среднего тригонометрической функцией (5/п) с учетом известного 
смещения фазы. 

Казалось бы, чего же ещё желать исследователю? Модельный график неплохо согласуется с 
данными наблюдений, выведен интегральный показатель активности (величина скользящего 
среднего)... Задача решена, не так ли? 

Оказывается - нет. В ходе исследования была совершена немалая методическая ошибка, а 
именно: в модель было привнесено извне желаемое толкование исследуемой динамики... то 
самое «скользящее среднее». Кроме того, часть важной информации была отброшена предпо¬ 
ложением, что «там нет никакой другой причины колебаний, кроме периодичности, задавае¬ 
мой тригонометрической функцией». А в результате - желаемое было выдано за действитель¬ 
ное. Истинная причина осталась за пределами понимания, а ведь поиск (и открытие) такой 
причины - основная, если не главная задача исследователя. 

3. Новация: поиск внутренних закономерностей 

Начнем сначала. В исследуемой динамике солнечной активности изначально содержится ин¬ 
формация, которую надо извлечь и, по возможности, наглядно продемонстрировать. Однако 
на этом пути возникает ряд сложностей. 

Проблема первая: для изучения внутренних закономерностей ряда чисел Вольфа не надо при¬ 
влекать (во всяком случае - пока, руководствуясь незыблемым принципом «бритвы Оккама») 
никаких других сведений. Значит, решение надо искать через возможные преобразования 
представления самого ряда. 























































































































Проблема вторая: для подобного исследования недостаточно простого, «одномерного» пред¬ 
ставления ряда... Ѵоіепз-поіепз, придется переходить от плоских графиков вида «годы» - «зна¬ 
чения» к объемным, «многомерным» фигурам. 

Проблема третья: до сих пор было непонятно, что это за фигуры, и как их сначала построить, а 
потом - и увидеть. 

Представим исходный ряд среднегодовых чисел Вольфа (ѵу) как последовательность ( іѴі 70 о, 
ѵѵіуоі, ..., іѵ 2 оі 9 )• Для того, чтобы учесть предысторию, заменим в этом ряду каждое одномерное 
значение щ вектором Щ = { щ , щ.і , ..., щ . 2 \}, то есть дополним значение числа Вольфа теку¬ 
щего года предшествующими (за 21 год) значениями. Видимо, не лишним будет напоминание 
о том, что «истинной» продолжительностью цикла солнечной активности следует считать не 11, 
а именно 22 года, с учетом изменения полярности магнитного поля Солнца 11 . 

В результате получим матрицу, строки которой (объекты, в принятых терминах многомерной 
статистики) будут соответствовать точкам в многомерном (размерности 22) пространстве. Вся 
совокупность таких точек образует опорные точки некой кривой, описывающей динамику ис¬ 
следуемого процесса. Осталось немного - «увидеть» эту кривую. Хотя, конечно, «видеть» в 22- 
мерном пространстве человеку (с его привычным 2,5-мерном представлением об окружающем 
мире) довольно-таки затруднительно. Однако заметим - с преобразованием исходного «плос¬ 
кого» ряда Вольфа к его многомерному аналогу первые две проблемы оказались преодолены. 

3.1. К НАГЛЯДНОСТИ ЧЕРЕЗ ПОНИЖЕНИЕ РАЗМЕРНОСТИ 

Очевидно, чтобы «увидеть» ход процесс в многомерном пространстве, надо предварительно 
преобразовать это пространство к такому виду, где количество значимых измерений сущест¬ 
венно меньше количества исходных признаков (столбцов матрицы данных). Для понижения 
размерности пространства существует несколько методов, из которых для наших целей вполне 
годится метод главных компонент 12 , фактически «изобретенный» Карлом Пирсоном 13 опирав¬ 
шимся на разработанный Джеймсом Сильвестером 14 математический аппарат. 

Суть метода применительно к нашей задачи: в исходной матрице данных информация «разма¬ 
зана» между признаками, образующими неортогональное пространство. Другими словами, ис¬ 
ходные признаки как-то коррелируют между собой, в результате чего вариации одного призна¬ 
ка оказываются связанными с вариациями других. При этом информативность каждого из при¬ 
знаков (определяемая через их дисперсию) практически не отличается от информативности 
любого из остальных признаков. 

Процедура преобразования («поворота») пространства методом главных компонент такова: 

1) вычисление матрицы парных корреляций между значениями признаков, 

2) разложение полученной матрицы корреляций на матрицу собственных векторов и век¬ 
тор собственных чисел, 

3) «поворот» исходной матрицы данных перемножением ее на матрицу собственных век¬ 
торов и вектор собственных чисел (упорядоченные по убыванию значений собственных 
чисел). 


11 Ю.И.Витинский, М.Копецкий, Г.В.Куклин «Статистика пятнообразовательной деятельности Солнца» — М: Наука, 
1986 

12 англ. РгіпсіраІ сотропепг апаіузіз, РСА 

13 К.Реагзоп «Оп Ііпез апб ріапез о^ сіозезг Яг го зузгетз ог роіпгз іп зрасе» - РЫІозорЫсаІ Ма§агіпе 2:559-572 (1901) 

14 .У.буіѵезгег «Оп гНе гесіисгіоп ог а Ьіііпеаг циапгіс ог гЬе пгЬ огсіег го гНе гогт ог а зит ог п ргосіисгз Ьу а сІоиЫе 
оггЬо§опаІ зиЬзгігигіоп» - Меззеп§ег ог МагИетагісз, 19 (1889), 42-46 



Таким образом, ранее существовавшие объекты получают в преобразованном пространстве 
новые координаты (соответствующие новым, вычисленным, осям) - но их относительное взаи¬ 
морасположение не меняется. 

Основным свойством пространства, описываемого матрицей собственных векторов, является 
его ортогональность; вариации по любому из собственных векторов (ортов) не влияют на ва¬ 
риации по остальным ортам. Дисперсия же распределяется между ортами неравномерно- 
«старшие» орты (имеющие меньшие номера) концентрируют больше дисперсии, чем «млад¬ 
шие». Распределение дисперсий по осям преобразованного пространства соответствует полу¬ 
ченным в ходе преобразования собственным числам, в чем легко убедиться, распечатав соот¬ 
ветствующие показатели входе вычислений. 

Выбирая из всего полученного набора главных компонент лишь те «старшие», суммарная дис¬ 
персия которых меньше заданного порогового значения, можно исключить из рассмотрения те 
«младшие», информативность которых не повлияет на дальнейшие исследования. 

3.2. Траектория процесса и её визуализация 

Как уже говорилось, переход от «плоского» ряда чисел Вольфа к его многомерному аналогу, 
опирающемуся на представления годичных точек как векторов Щ = { щ , нѵі ,..., Щ-п }, создает 
вместо простого графика «годы» - «числа» многомерную кривую, опирающуюся на годичные 
точки. Для визуализации этой кривой достаточно выполнить ее проецирование на фазовые 
плоскости (главные компоненты), образованные ортами - собственными векторами. Важно от¬ 
метить, что никаких новых предположений о характере изучаемого процесса при таком преоб¬ 
разовании не делалось, и никакие «желаемые» зависимости в исходные данные не вносились. 

Для получения первого (наиболее очевидного) представления о динамическом процессе сол¬ 
нечной активности достаточно «нарисовать» хотя бы три проекции, на плоскости І-ІІ, ІІ-ІІІ и І-ІІІ 
главных компонент. 

Отдельно следует рассмотреть «вклады» годов (строго говоря, лагов - «запаздываний») в полу¬ 
ченные собственные векторы. Сумма квадратов всех вкладов для любого собственного вектора 
равна единице. Рассматривая вклады лагов в порядке убывания их абсолютного значения, 
можно оценить накопленную сумму; достижение заданного порога значимости можно воспри¬ 
нимать как сигнал к прекращению дальнейшего рассмотрения вкладов, и при толковании при¬ 
роды вариаций по данной оси преобразованного пространства ограничиться только набором 
уже отмеченных вкладов. 

По всей видимости, проекция исходных признаков на фазовые плоскости главных компонент 
также имеет смысл, поскольку отражает (в данном случае) тенденции многолетних, повторяю¬ 
щихся изменений в характере влияния факторов изучаемого процесса на поведение динамиче¬ 
ской системы. 

В изучение динамических процессов (временных рядов) методом главных компонент сущест¬ 
венный вклад внес Вадим Ефимов 15 , первым предложивший 16 , в частности, методику последо¬ 
вательного соединения точек-объектов на дискретных проекциях непрерывной кривой. Автор 
выражает В.М.Ефимову глубокую искреннюю благодарность за возможность приобщения к 
идеям и методам многомерной статистики. 


15 В.М.Ефимов «Обработка временных рядов методом главных компонент» - науч.-техн.бюл./ВАСХНИЛ, Сиб.отд- 
ние, вып.22, с.32-40. - Новосибирск: 1984 

16 Ю.К.Галактионов, частное сообщение (1985) 



4. Результаты 


Вычисления проводились по программе \л/оІ/-поѵа.іас (см. Приложение П2.3). Полученные в 
ходе расчетов значения собственных чисел демонстрируют радикальное снижение размерно¬ 
сти при переходе от исходного пространства признаков к пространству главных компонент 

(Табл.З): 


Табл.З. Собственные числа [дисперсии компонент) для 22-летнего лага 



I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

XI 

А 

7,25 

7,03 

2,88 

1,23 

0,97 

0,91 

0,39 

0,32 

0,32 

0,17 

0,11 

% 

32,97 

31,95 

13,07 

5,60 

4,42 

4,15 

1,79 

1,46 

1,45 

0,77 

0,50 

2% 

32,97 

64,92 

77,99 

83,59 

88,02 

92,16 

93,96 

95,42 

96,87 

97,64 

98,15 


В приведенной таблице столбцы соответствуют главным компонентам (приведены 11 из 22, ос¬ 
тальные компоненты не имеют существенного значения), Л - значения собственных чисел, % - 
доля дисперсии, относящейся к соответствующей главной компоненте (в процентах), Е%- на¬ 
копленная доля дисперсии (в процентах). Жирным шрифтом выделены показатели накоплен¬ 
ной доли дисперсии для І-ѴІ главных компонент, в сумме описывающих не менее 90% общей 
вариации (порог информативности 0,90). 

Очевидно, что основная доля дисперсии (соответственно, значительная доля информации о 
процессе) приходится на І-ІІІ главные компоненты. Посмотрим, как ведет себя процесс в новых, 
ортогональных координатах. 



Рис.З. Проекция объектов (годы) на плоскость І-ІІ главных компонент 


На проекции годы отмечены точками, цветом выделены два периода, представляющие опре¬ 
деленный исторический (и не только) интерес. Голубой линией отмечены годы, относящиеся к 





































так называемому «минимуму Далтона», во многом совпадающим с годами Великой Француз¬ 
ской революции и последующими потрясениями, охватившими значительную часть Европы. 
Красной линией выделены годы после 2000-го, события которых известны большинству читате¬ 
лей не понаслышке. 

Налицо явная цикличность процесса. Но только ли об этой, «лежащей на поверхности» циклич¬ 
ности говорят нам числа Вольфа? Взглянем на тот же процесс - «вид сбоку». 



Рис.4. Проекция объектов на плоскость І-ІІІ (слева) и П-Ш (справа) главных компонент 


Траектория выглядит так, словно воображаемая точка, характеризующая процесс во времени, 
«пробегает» по гипотетической «воронке», причём заметно - чем выше точка располагается по 
вертикальной координате (она же - III главная компонента), тем больше диаметр описываемой 
ею окружности. Что это, если не обобщённый показатель энергетики процесса? 

Проверим, насколько смещения по вертикали соответствуют радиусу описываемой окружно¬ 
сти. Как известно, радиус окружности /? = (х 2 +у 2 ) 1/2 , где хи у- координаты точки в плоскости (I и 
II главных компонент). Вычислим этот радиус для каждой точки, и определим корреляцию ме¬ 
жду ним и соответствующим значением III главной компоненты. А заодно проверим, насколько 
этот показатель соответствует уже использовавшемуся значению «скользящего среднего» (про¬ 
грамма шоЦ-соп.іас, см. Приложение П2.4). 

Корреляция между вычисленными значениями радиуса и значениями III компоненты составля¬ 
ет 0,375, в то время, как корреляция между значениями скользящего среднего и значениями 
радиуса и III компоненты составляет, соответственно, 0,018 и 0,073. Очевидно, что применение 
«метод скользящего среднего» для исследования внутренних закономерностей динамики сол¬ 
нечной активности лишено всякого смысла. 

С другой стороны, обнаруживаются значительные корреляции между радиусом Кщ, вычислен¬ 
ным по значениям I и II главных компонент, и радиусом Я| Ѵ ,ѵ,ѵі, вычисленным по значениям 
компонент IV, V и VI (надо отметить, что в соответствии с данными Табл.З эти компоненты яв¬ 
ляются значимыми для статистического описания исследуемого процесса). Значение коэффи¬ 
циента корреляции между этими параметрами составляет 0, 465, что уже превышает корреля¬ 
цию между значениями радиуса /?і,ц и значениями III компоненты. Кстати, коэффициент корре¬ 
ляции между значениями III компоненты и радиусом /?іѵ,ѵ,ѵі составляет 0,484. Это неоспоримо 
свидетельствует о наличии существенных внутренних зависимостей в динамике солнечной ак¬ 
тивности, выражаемой даже таким косвенно-опосредованным показателем, как числа Вольфа. 

Следовательно, в ходе проводимого анализа данных необходимо уделить внимание поведе¬ 
нию траектории исследуемого процесса в фазовом пространстве IV, V и VI компонент. 













































Рис.5. Проекция объектов на плоскость: ІѴ-Ѵ, ІѴ-ѴІ и Ѵ-ѴІ главных компонент 

Проявления периодических процессов заметны, но дать им какое-то однозначное толкование 
на основании проекций (Рис.5) невозможно. При этом вполне отчетливо выделяются годы «ми¬ 
нимума Далтона», что должно навести на мысль о его неслучайной природе. Удовлетворимся 
тем, что сильно коррелирующую с общей амплитудой процесса комбинацию IV, V и VI компо¬ 
нент обозначим как «фактор X», и посмотрим на поведение его и энергетического «факто¬ 
ра Е» на фоне уже знакомого ряда среднегодовых чисел Вольфа (Рис.б). 

300 т- 



Черным цветом показан график изменения среднегодовых чисел Вольфа, фиолетовым - изме¬ 
нения «фактора Е», красным - изменения «фактора X». Ординаты для графиков факторов ус¬ 
ловны, их значения выражают только относительные колебания. 

Небезынтересно будет взглянуть и на поведение исходных признаков (лагов лет) в построен¬ 
ном методом главных компонент фазовом пространстве. В типичном случае (лаг 22 года) гра¬ 
фики процесса получаются достаточно «короткими», поэтому для начала увеличим лаг до 69 
лет (программа \л/оІ/-сІеерфас, см. Приложение П2.5). 

Подтверждением допустимости использования такого лага лет служат значения «старших» соб¬ 
ственных чисел, причем не сами значения, а накопленная доля общей дисперсии первых шести 
компонент. Как и в предыдущих расчетах, значимыми остаются компоненты І-ѴІ. Воспользуем¬ 
ся распечаткой матрицы собственных векторов, и используем её в качестве основы для по¬ 
строения проекций признаков (лагов лет) на фазовые плоскости главных компонент. 

























































































































Первые три проекции показывают поведение лагов лет в пространстве І-ІІІ главных компонент, 
и отражают уже привычную цикличность, сопряженную с переменной общей энергетикой ди¬ 
намического процесса солнечной активности (Рис.7): 



Рис.7. Фазовые портреты лага лет в пространстве І-ІІ, І-ІІІ и ІІ-ІІІ главных компонент 


Конечно, после проведенного исследования поведения траектории объектов-лет в пространст¬ 
ве первых трёх главных компонент такие фазовые портреты оказываются вполне ожидаемыми. 
Уже понятно, что процесс имеет ярко выраженную цикличность. С другой стороны, эта циклич¬ 
ность прочно связана с изменением энергетической составляющей - и вот здесь фазовые порт¬ 
реты лагов лет в пространстве III компоненты наводят на мысль, что такие изменения вовсе не 
случайны, что существует вполне устойчивый механизм, систематически «разворачивающий» 
плоскость вращения лагов лет относительно III главной компоненты. 

Посмотрим ещё глубже, не отыщем ли каких-нибудь доселе скрытых закономерностей и там... 



Рис.8. Фазовые портреты лага лет в пространстве ІѴ-Ѵ, ІѴ-ѴІ и Ѵ-ѴІ главных компонент 


Выходит, кроме уже известной «воронки», связанной с общей динамикой появления солнеч¬ 
ных пятен, есть и ещё... И не одна, а сразу две. Воображаемая точка лага лет описывает в про¬ 
странстве ІѴ-ѴІ главных компонент траекторию, отчасти напоминающую две «воронки», соеди¬ 
ненные узкими частями через «перемычку». Соответственно, объединяющий эффекты ІѴ-ѴІ 
компонент «фактор X» сначала как бы «сворачивается» до минимальных амплитуд по одной из 
связанных «воронок», потом проходит без каких-либо резких изменений через «перемычку», а 
потом «раскручивается» по второй «воронке»... И так снова и снова, в течение многих лет. 

Выбрав в качестве лага 69 лет, видим уже как минимум четырехкратное воспроизведение этого 
сложного, глубоко скрытого квазициклического процесса за всю историю наблюдений солнеч¬ 
ной активности и фиксации среднегодовых значений чисел Вольфа. Так что это - не случай¬ 
ность, и не следствие «ошибок метода». 

Это существует на самом деле, и - надеюсь! - в недалеком будущем прояснится и физическая 
природа продемонстрированной закономерности. 




























5. Возражения 


Конечно, возражения будут. Слишком долго на динамику солнечной активности смотрели по¬ 
верхностно, ограничиваясь построением таких простейших моделей, которые, как бы воспро¬ 
изводя внешние эффекты изменений в наблюдаемых проявлениях активности, позволяли 
удовлетвориться тривиальными объяснениями. Ну, есть цикличность в появлении и исчезнове¬ 
нии пятен на Солнце. Ну, посчитали продолжительность цикла - 11 лет (примерно). Попытались 
как-то связать 17 эту цикличность с её земными проявлениями. На том, по сути, и успокоились. 

Конечно, с опубликованием этой работы претензии автору будут предъявляться как по поводу 
метода выявления закономерностей в проявлениях солнечной активности, так и по части избы¬ 
точности описания обнаруженных явлений. 

Рассмотрим влияние методики учета предыстории текущего состояния процесса. Немного из¬ 
меним программу вычислений, заменив реальные значения среднегодовых чисел Вольфа слу¬ 
чайными значениями, подчиняющимися тем же законам распределения (сгенерируем их с те¬ 
ми же параметрами, что и у чисел реального ряда, и/ тіп =0, и/ тах =269,2), сохранив длину ряда. 
Вычисления произведены программой и/о//-гопс/.уос (см. Приложение П2.6). 

Прежде всего обратим внимание на полученные собственные числа (Табл.4). Налицо более 
«плавное» уменьшение информационной значимости компонент, причем концентрации ин¬ 
формации в «старших» компонентах не происходит. 


Табл.4. Собственные числа (дисперсии компонент) для случайного ряда 



I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

XI 

А 

1,49 

1,46 

1,43 

1,18 

1,11 

1,09 

1,09 

1,08 

1,06 

1,00 

0,96 

% 

6,77 

6,64 

6,48 

5,35 

5,04 

4,97 

4,95 

4,91 

4,83 

4,55 

4,37 

і% 

6,77 

13,40 

19,88 

25,23 

30,28 

35,25 

40,20 

45,11 

49,94 

54,49 

58,86 


Поведение объектов («лет») еще более характерно (Рис.9). Нигде не наблюдается той циклич¬ 
ности, которая была наглядно продемонстрирована ранее. 



Рис.9. Проекции объектов случайного ряда на плоскости І-ІІ, І-ІІІ и ІІ-ІІІ главных компонент 


Очевидно, что подобная картина будет наблюдаться и во всех следующих проекциях... 

Настала пора сравнить две модели - традиционную (метод скользящего среднего) и новацион¬ 
ную (метод главных компонент). Сторонники традиционной модели вправе заявить: наша мо¬ 
дель учитывает и краткосрочную, и долгосрочную динамику солнечной активности. Хорошо, 
попробуем использовать данные, полученные традиционным моделированием, как исходные 
для анализа методом главных компонент. 
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А.Л.Чижевский «Земное эхо солнечных бурь» — М: Мысль, 1973 

































Для вычислений воспользуемся программой и /оЦ-сотр.]ас (см. Приложение П2.7). Сначала по¬ 
лучим модельные значения традиционным способом, а потом применим к ним стандартную 
схему вычисления главных компонент с привычным 22-летним лагом. 

Полученные дисперсии компонент (Табл.5) показывают, что полноценно-значимыми являются 
лишь компоненты І-ІІІ. Нно всё же бросим взгляд на поведение объектов-лет в пространстве до 
VI главной компоненты (в том подмножестве, значимость которого достигает 99%). 

Табл.5. Собственные числа (дисперсии компонент) для традиционной модели 



I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

XI 

А 

9,04 

9,00 

2,39 

0,78 

0,34 

0,29 

0,09 

0,03 

0,01 

0,01 

0,01 

% 

41,09 

40,90 

10,87 

3,56 

1,52 

1,31 

0,42 

0,14 

0,06 

0,05 

0,03 

2% 

41,09 

81,99 

92,87 

96,43 

97,96 

99,26 

99,69 

99,83 

99,89 

99,94 

99,97 


И вот, наконец - фазовые портреты поведения траектории объектов (лет) в пространстве наи¬ 
более значимых компонент (Рис.10): 



Рис.10. Фазовые портреты традиционной модели в плоскости І-ІІ, І-ІІІ и ІІ-ІІІ главных компонент 


Похоже на рисунки, относящиеся к новационной модели?-да, конечно. Но какое-то оно всё- 
слишком гладкое, слишком геометрическое. Да и следующие компоненты какие-то малоин¬ 
формативные (Рис.11): 



Рис.11. Фазовые портреты традиционной модели в плоскости ІѴ-Ѵ, ІѴ-ѴІ и Ѵ-ѴІ главных компонент 


Применение принудительной фильтрации в виде метода скользящего среднего привело к утра¬ 
те представлений о зависимостях второго порядка. Жаль, но использование традиционной мо¬ 
дели делает исследователя подслеповатым, и при всём желании он не разглядит эффектов, 
лежащих глубже простых, тривиальных пульсаций, описываемых тригонометрическими функ¬ 
циями вкупе со «скользящим средним». При таком традиционном подходе не стоит даже и 
заикаться о том, чтобы посмотреть на поведение лагов лет в пространстве главных компонент, 
поскольку из модельных данных эти эффекты исчезли. 




















































Безусловно, множество возможных возражений не ограничивается двумя, рассмотренными 
выше. Однако представляется справедливым попросить (на будущее) всех критиков предла¬ 
гаемого подхода заранее, до предъявления автору претензий, самостоятельно выполнить вы¬ 
числительную проверку своих аргументов, благо весь инструментарий для этого доступен, от¬ 
крыт и (какой ужас!) совершенно бесплатен. 

6. Заключение 

Итак, для детального исследования были взяты данные, известные человечеству на протяже¬ 
нии сотен лет. К этим данным был применен метод, известный (как подход и математический 
аппарат) уже более ста лет. Единственной новацией стало применение инструмента, позво¬ 
ляющего исследователю не сосредотачиваться на тонкостях машинной реализации тех или 
иных вычислений, но направлять свой интеллект и любознательность на изучение самого явле¬ 
ния, благо создание и выполнение схем обработки данных (ЯКОБИ-программ) с помощью 
калькулятора многомерной статистики не требует сверхъестественных талантов и усилий. 

Обнаруженные зависимости еще ждут своего объяснения, еще только предстоит найти и сфор¬ 
мулировать физическую причину таких проявлений динамики солнечной активности. Скорее 
всего, для глубокого изучения явлений понадобится провести аналогичную обработку более 
детальной информации (среднемесячных чисел Вольфа, например). Неизбежно привлечение 
данных объективных измерений - однако для таких параметров период наблюдений окажется 
существенно меньшим по продолжительности. 

По всей видимости, откроется влияние выделенных в данном исследовании «фактора Е» и 
«фактора X» не только на физические процессы, но и на биосферу Земли вообще, и историю 
человечества в частности (особенно интересным видится сопоставление лет, «отмеченных» вы¬ 
сокими значениями «фактора X» с рождением, жизнью и творчеством ярких личностей). 

Разумеется, данное исследование не претендует на определение «истины в последней инстан¬ 
ции», в отечественной традиции существует не один способ усомниться 18 в правомочности по¬ 
добных " штудий". Тем не менее, описываемые закономерности существуют и, как бы кому-то 
не хотелось избегнуть их упоминания (тем более - называния настоящим, подлинным име¬ 
нем), с этим придется смириться. Те же благодарные читатели, которые рискнут продолжить 
начатое этой работой изучение доселе скрытых от внимания общества явлений, отныне имеют 
в своем распоряжении все необходимые данные и инструменты для их анализа. Для бескоры¬ 
стного добывания знаний об окружающем нас Мире. 

Ради всего человечества 19 . 


18 В.Семи-Булатов «Письмо к ученому соседу» (цит. по: А.П.Чехов - Поли. собр. соч. и писем в 30-ти томах) - 
Том 1, с.11-17. - М: «Наука», 1974 

19 англ. «Рог аІІ тапкіпсі» - слова, произнесённые Нилом Олденом Армстронгом (англ. ЫеіІ АІсІеп Агпл5І:гоп§; 5 ав¬ 
густа 1930 - 25 августа 2012), первым человеком, шагнувшим на поверхность Луны, 20 июля 1969 года (автор 
был косвенным свидетелем этого события, слушая прямой эфир «Голоса Америки» недалеко от Новосибирска) 



Приложения 


П1. Среднегодовые числа Вольфа (1700-2019 годы) 

Приводимая таблица содержит значения 20 среднегодовых чисел Вольфа, для компактности 
представления сгруппированные по десятилетиям; в заголовке указано смещение года от на¬ 
чала соответствующего десятилетнего периода. 


Табл.ПІ.1. Среднегодовые числа Вольфа по десятилетиям за период 1700-2019 годов 


Год 

+0 

+1 

+2 

+3 

+4 

+5 

+6 

+7 

+8 

+9 

1700 

8,3 

18,3 

26,7 

38,3 

60,0 

96,7 

48,3 

33,3 

16,7 

13,3 

1710 

5,0 

0,0 

0,0 

3,3 

18,3 

45,0 

78,3 

105,0 

100,0 

65,0 

1720 

46,7 

43,3 

36,7 

18,3 

35,0 

66,7 

130,0 

203,3 

171,7 

121,7 

1730 

78,3 

58,3 

18,3 

8,3 

26,7 

56,7 

116,7 

135,0 

185,0 

168,3 

1740 

121,7 

66,7 

33,3 

26,7 

8,3 

18,3 

36,7 

66,7 

100,0 

134,8 

1750 

139,0 

79,5 

79,7 

51,2 

20,3 

16,0 

17,0 

54,0 

79,3 

90,0 

1760 

104,8 

143,2 

102,0 

75,2 

60,7 

34,8 

19,0 

63,0 

116,3 

176,8 

1770 

168,0 

136,0 

110,8 

58,0 

51,0 

11,7 

33,0 

154,2 

257,2 

209,8 

1780 

141,3 

113,5 

64,2 

38,0 

17,0 

40,2 

138,2 

220,0 

218,2 

196,8 

1790 

149,8 

111,0 

100,0 

78,2 

68,3 

35,5 

26,7 

10,7 

6,8 

11,3 

1800 

24,2 

56,7 

75,0 

71,8 

79,2 

70,3 

46,8 

16,8 

13,5 

4,2 

1810 

0,0 

2,3 

8,3 

20,3 

23,2 

59,0 

76,3 

68,3 

52,9 

38,5 

1820 

24,2 

9,2 

6,3 

2,2 

11,4 

28,2 

59,9 

83,0 

108,5 

115,2 

1830 

117,4 

80,8 

44,3 

13,4 

19,5 

85,8 

192,7 

227,3 

168,7 

143,0 

1840 

105,5 

63,3 

40,3 

18,1 

25,1 

65,8 

102,7 

166,3 

208,3 

182,5 

1850 

126,3 

122,0 

102,7 

74,1 

39,0 

12,7 

8,2 

43,4 

104,4 

178,3 

1860 

182,2 

146,6 

112,1 

83,5 

89,2 

57,8 

30,7 

13,9 

62,8 

123,6 

1870 

232,0 

185,3 

169,2 

110,1 

74,5 

28,3 

18,9 

20,7 

5,7 

10,0 

1880 

53,7 

90,5 

99,0 

106,1 

105,8 

86,3 

42,4 

21,8 

11,2 

10,4 

1890 

11,8 

59,5 

121,7 

142,0 

130,0 

106,6 

69,4 

43,8 

44,4 

20,2 

1900 

15,7 

4,6 

8,5 

40,8 

70,1 

105,5 

90,1 

102,8 

80,9 

73,2 

1910 

30,9 

9,5 

6,0 

2,4 

16,1 

79,0 

95,0 

173,6 

134,6 

105,7 

1920 

62,7 

43,5 

23,7 

9,7 

27,9 

74,0 

106,5 

114,7 

129,7 

108,2 

1930 

59,4 

35,1 

18,6 

9,2 

14,6 

60,2 

132,8 

190,6 

182,6 

148,0 

1940 

113,0 

79,2 

50,8 

27,1 

16,1 

55,3 

154,3 

214,7 

193,0 

190,7 

1950 

118,9 

98,3 

45,0 

20,1 

6,6 

54,2 

200,7 

269,2 

261,7 

225,1 

1960 

159,0 

76,4 

53,4 

39,9 

15,0 

22,0 

66,8 

132,8 

150,0 

149,3 

1970 

148,0 

94,4 

97,6 

54,1 

49,2 

22,5 

18,4 

39,3 

131,0 

220,1 

1980 

218,8 

198,8 

162,3 

91,0 

60,5 

20,6 

14,8 

33,9 

123,0 

211,1 

1990 

191,8 

203,3 

133,0 

76,1 

44,9 

25,1 

11,6 

28,9 

88,3 

136,3 

2000 

173,8 

170,3 

163,6 

99,3 

65,3 

45,8 

24,7 

12,6 

4,2 

4,8 

2010 

24,9 

80,8 

84,5 

94,0 

113,3 

69,8 

39,8 

21,7 

7,0 

3,6 


20 В работе использованы значения среднегодовых чисел Вольфа, опубликованные на веб-сайте Мирового Центра 
Данных (511-50 ба1:а/іта§е, КоуаІ ОЬзегѵаІогу оі ВеІ^ішп, Вгиззеіз), Ы:1:р://5ІсІс.Ье/5ІІ5о/сІа1=аГіІе5 




















































П2. ЯКОБИ-программы обработки данных 

Приводится только содержательная часть программ; директивы управления протоколом, про¬ 
межуточной печати и проецирования объектов/признаков на фазовые плоскости не показаны. 

П2.1. Программа ѵѵо1Моп§.)ас - расчет корреляций между текущим и предшествующими 
годами наблюдений солнечной активности за длительный период, исходные дан¬ 
ные ѵѵоИ'.зісі - ряд среднегодовых чисел Вольфа 

* читать данные * мэі-р 

* ввести объекты с номерами * ][ ][ 

* создать признаки * [1-149] 

* копировать значения признака * ][ [0] ]0[ 

* сдвинуть столбцы * ][ ]0[ +1 

* исключить объекты с отсутствием данных * ][ ][ 

* центрировать столбцы * ][ ][ 

* нормировать столбцы * ][ ][ 

* вычислить главные компоненты * ][ ][ 1 

* печатать ковариации(корреляции) * [0] [0-9] 

* * [0] [10-19] 

* * [0] [20-29] 

* * [0] [30-39] 

* * [0] [40-49] 

* * [0] [50-59] 

* * [0] [60-69] 

* * [0] [70-79] 

* * [0] [80-89] 

* * [0] [90-99] 

* * [0] [100-109] 

* * [0] [110-119] 

* * [0] [120-129] 

* * [0] [130-139] 

* * [0] [140-149] 

* конец * 

П2.2. Программа ѵѵоІГ-Ігай.іас - модельный расчет динамики солнечной активности на 
основе тригонометрической функции и скользящего среднего, исходные данные 
ѵѵоІГ.зМ - ряд среднегодовых чисел Вольфа 

* читать данные * моі-р 

* ввести объекты с номерами * ][ ][ 

* создать признаки * [1-10] 

* копировать значения признака * ][ [0] ]0[ 

* сдвинуть столбцы * ][ [1-10] +1 

* создать признаки * [77,88,99] 

* присвоить значения * ][ [77] 0 

* изменить значения с шагом * ]1700[ [77] +0.5754 

* сдвинуть столбцы * ][ [77] +3 

* стандартная функция * ][ [77] 5іп 

* действие с константой * ][ [77] +1 

* действие с признаком * ][ [88] : [0-10] 

* присвоить значения * ][ [99] 1 

* действие с признаком * ][ [99] * [77,88] 

* записать объекты с номерами * моі-р-іігасі ][ [0,88,99] 

* конец * 



П2.3. Программа лѵоІГ-поѵа^ас - расчет главных компонент с учетом 22-летнего лага, ис¬ 
ходные данные ѵѵоі^.зігі - ряд среднегодовых чисел Вольфа 

* читать данные * моі-р 

* ввести объекты с номерами * ][ ][ 

* создать признаки * [1-21] 

* копировать значения признака * ][ [0] ]0[ 

* сдвинуть столбцы * ][ ]0[ +1 

* исключить объекты с отсутствием данных * ][ ][ 

* центрировать столбцы * ][ ][ 

* нормировать столбцы * ][ ][ 

* вычислить главные компоненты * ][ ][ 1 

* печатать дисперсии компонент * ][ 0.90 

* записать объекты с номерами * іл/оі-р-поѵа ][ [1-6] 

* конец * 

П2.4. Программа шоИ-согг.)ас - расчет корреляций факторов новационной модели и 
скользящего среднего традиционной модели, исходные данные ѵѵоІГ.зЫ- ряд 
среднегодовых чисел Вольфа; как подпрограммы вызываются ранее представлен¬ 
ные программы ѵѵоІМгагі.іас и шоіі'-поѵа.іас 

* выполнить * моІ-Р-'ЬгасІ 

* записывать данные * іл/оі-р-аѵег 

* записать объекты с номерами * ][ [0,88,99] 

* закончить запись * 

* очистить память данных * 

* выполнить * моі-р-поѵа 

* исключить фрагмент матрицы данных * [] ]1-б[ 

* создать признаки * [12,46] 

* присвоить значения * ][ [12] 0 

* стандартная функция * ][ [1-2] Бцг 

* действие с признаком * ][ [12] + [1-2] 

* стандартная функция * ][ [12] Бцгі: 

* присвоить значения * ][ [46] 0 

* стандартная функция * ][ [4-6] Бцг 

* действие с признаком * ][ [46] + [4-6] 

* стандартная функция * ][ [46] Бцгі: 

* записывать данные * моі-р-сотр 

* записать объекты с номерами * ][ ][ 

* закончить запись * 

* очистить память данных * 

* создать признаки * [0,1-6,12,46,88,99] 

* читать данные * іл/оі-р-аѵег 

* ввести объекты с номерами * ][ [0,88,99] 

* читать данные * іл/оі-р-сотр 

* ввести объекты с номерами * ][ [1- 6,12,46] 

* записать объекты с номерами * іл/оі-р-зитт ][ ][ 

* исключить фрагмент матрицы данных * [] ]0,3,12,46,88[ 

* исключить объекты с отсутствием данных * ][ ][ 

* центрировать столбцы * ][ ][ 

* нормировать столбцы * ][ ][ 

* вычислить главные компоненты * ][ ][ 

* печатать ковариации(корреляции) * ][ ][ 

* конец * 



П2.5. Программа ѵѵоІГ-сіеер.іас - вычисление собственных векторов при заведомо боль¬ 
шом лаге лет, исходные данные ѵѵо1Г.5Ігі - ряд среднегодовых чисел Вольфа 

* читать данные * моі-р 

* ввести объекты с номерами * ][ [0] 

* создать признаки * [1-68] 

* копировать значения признака * ][ [0] ]0[ 

* сдвинуть столбцы * ][ ]0[ +1 

* исключить объекты с отсутствием данных * ][ ][ 

* центрировать столбцы * ][ ][ 

* нормировать столбцы * ][ ][ 

* вычислить главные компоненты * ][ ][ 

* печатать дисперсии компонент * ][ 0.7 

* печатать собственные векторы * ][ ][ 

* конец * 

П2.6. Программа шоН-гапс1.)ас - расчет главных компонент с учетом 22-летнего лага, ис¬ 
ходные данные - ряд случайных чисел с теми же параметрами, что и исходный ряд 
среднегодовых чисел Вольфа 

* читать данные * мэі-р 

* ввести объекты с номерами * ][ ][ 

* присвоить случайные значения * ][ [0] К 0 269.2 

* создать признаки * [1-21] 

* копировать значения признака * ][ [0] ]0[ 

* сдвинуть столбцы * ][ ]0[ +1 

* исключить объекты с отсутствием данных * ][ ][ 

* центрировать столбцы * ][ ][ 

* нормировать столбцы * ][ ][ 

* вычислить главные компоненты * ][ ][ 

* печатать дисперсии компонент * ][ 0.75 

* записать объекты с номерами * іл/оі-р-гапсі ][ [1-6] 

* конец * 

П2.7. Программа ѵѵоіі-сотр.іас - расчет главных компонент с учетом 22-летнего лага, 
исходные данные - модельный ряд, полученный из исходного ряда среднегодовых 
чисел Вольфа 

* выполнить * моІ-Р-'ЬгасІ 

* копировать значения признака * ][ [99] [0] 

* исключить фрагмент матрицы данных * [] ]0[ 

* создать признаки * [1-21] 

* копировать значения признака * ][ [0] ]0[ 

* сдвинуть столбцы * ][ ]0[ +1 

* исключить объекты с отсутствием данных * ][ ][ 

* центрировать столбцы * ][ ][ 

* нормировать столбцы * ][ ][ 

* вычислить главные компоненты * ][ ][ 

* печатать дисперсии компонент * ][ 0.9 

* записать объекты с номерами * іл/оі-р-сотр ][ [1-6] 

* конец * 



